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Abstract 

We investigate the crumpling transition on crystalline random surfaces 
with extrinsic curvature on lattices up to 64 2 . Our data are consistent 
with a second order phase transition and we find correlation length critical 
exponent v = 0.89 ± 0.07. The specific heat exponent, a = 0.2 ± 0.15, 
is in much better agreement with hyperscaling than hitherto. The long 
distance behaviour of tangent-tangent correlation functions confirms that 
the so-called Hausdorff dimension is djj = oo throughout the crumpled 
phase. 



The crystalline random surface plays an important role in a number of areas. 
Besides its obvious condensed matter applications Jl|, ||, it is also of relevance in 
investigations of the non-perturbative regularization of strings and the interaction 
of quantum gravity with matter in two dimensions 0. The model, hereafter 
referred to as "model 1" , is defined on a triangulated surface which forms a two- 
dimensional lattice of fixed topology; the lattice sites i have coordinates X(i) in 
the three-dimensional embedding space and the action is 

s = \ E ( x « - x 0')) 2 - * E ft A.ftA' (i) 

1 <ij> A A' 

where < ij > denotes the link fron i to j and tiA.A' are the unit normal vectors of 
the triangles on either side of the link. The second term in S, called the extrinsic 
curvature, tends to make the surface smoother. Expectation values of operators 
are defined by 

(•) = ^n/ rf3x (^)^(E X (j))(-)e- 5 (2) 

1 i j 

where Z\ = (1) and the delta function suppresses the translational zero mode. 
When k = the model is trivial and the surface has mean square extent 

^ = ^E X « 2 -logiV (3) 

i 

Introducing the so-called Hausdorff dimension, dn, through B? ~ N 2 / dH the 
model at K = has du = oo; this means that the surface is all crumpled up. 
There has been amassed a large body of evidence @, f|, || that at large enough k 
the surface becomes smooth with du = 2 and that at k = n c there is a phase tran- 
sition, called the crumpling transition, which separates these two regimes. The 
evidence suggests that the phase transition is second order but different mea- 
surements of the critical exponents ||, |[] have yielded different results, although 
usually consistent within the rather large estimated errors. Studies of strip ge- 
ometries have also been made to determine the central charge c at the transition 
[Bl 0. These strongly suggest that c e // < 1 but, because of the possibility of 
non-unitary behaviour, the interpretation of this is obscure. 

In the context of the string, the inclusion of an extrinsic curvature term was 
first proposed in [[|. It is believed that the discretized equivalent of Polyakov's 
rigid string, hereafter referred to as "model 2" , has partition function given by 

^ 2 = E /CO II / ^ 3x (^ 3 (E Mj))e~ S (4) 

T i j 

where J2t denotes the sum over all triangulations, T, of the surface and /(T) 
depends only on the triangulation. Summing over triangulations ||, |9], [K| is 
held to be equivalent to integrating over the metric in the continuum formulation 
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and it has been shown (TT| that the extrinsic curvature term in the action is 
necessary if there is to be scaling of the string tension in the continuum limit. 
Perturbation theory calculations in the continuum model do not predict a phase 
transition. However, there is substantial numerical evidence O, Il3l O, [TBI |l6l 



for some kind of transition in the lattice model. The most recent work [15|, [16 
finds weaker divergence in the specific heat than earlier simulations. 

According to conventional ideas about these models, the effective field theory, 
F 2 , governing the behaviour of model 2 at its critical point (should there be 
one) is the same theory as the effective field theory, F\ (governing the critical 
point of model 1), interacting with quantum gravity. In particular, should F\ be 
unitary, polynomial and have central charge c < 1 then the properties of F 2 and 
Fi should be related by the KPZ result |17|, |Tj|. However it is not clear either 



that Fi is unitary or polynomial or that it has c < 1; indeed one might hope 
that these models, which naively have more than one bosonic degree of freedom, 
provide a way of getting a c > 1 field theory interacting with gravity in a non- 
pathological manner and it is certainly the case that the naive continuum limit of 
(1) contains higher derivative terms and is not unitary while (1) itself is certainly 
not polynomial in the X fields. 

To settle the questions discussed above it is necessary to have high quality 
information about both models. We have been running a new simulation of model 
1 based on the Fourier accelerated Langevin algorithm used successfully in 0. 



The algorithm will be discussed in detail elsewhere [|19| but the renormalization of 
k (an inevitable consequence of discretizing the Langevin time) is different from 
that in || so that in the present work the critical value k c is not the same as in 
[Ej. Our aim has been to obtain much higher statistics than in [EJ but on similar 
lattice sizes; we have obtained data on lattices with toroidal boundary conditions 
and sizes of N = 16 2 , 16 x 32, 32 2 , 32 x 64, 64 2 . Using a basic Langevin step size 
of .002, on the largest lattice we have accumulated of order 2.5 x 10 6 Langevin 
updates at each k value studied, an improvement by a factor of O(10) over ||. 
Running will continue with larger sizes but the results so far are an advance over 
previous work and we present them here. 



As discussed at length in [EUL the tangent-tangent correlation functions 



yield a great deal of information about the system. Let 

t 1 (£i,6)=X(£i + l,6)-X(£i,&) (5) 

and define the correlation function G\\ by 

G|l(n) = (ti(£,0).ti(0,0)) (6). 

It is easy to show || that must become negative at large £ so that Grii(£o) = 

for some £o- Supposing the model has only one relevant length scale, the mass 
gap in the crumpled phase, m, satisfies 

m oc f^ 1 (7) 
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on an infinite size system. However, on a finite size system with L lattice spac- 
ings in the 1 direction the value of £o cannot exceed L/4 || and when the £o 
appropriate to L — > oo approaches L/4 we expect to see finite size effects. Fig.l 
shows £o extracted from G\\ on 16 2 , 32 2 and 64 2 systems and clearly illustrates 
the existence of these finite size effects. To unravel them requires some model; 
previous work || showed that the correlation functions in the crumpled phase 
were well described at long distances by supposing that the system is essentially 
a free field theory and that in momentum space 

< x ®- x <-^> « mmr^) (8) 

where L(p) is the lattice Laplacian. Making this assumption we can remove the 
finite size effects by fitting (8) to the data in the region where G\\ < Q 

The fit is discussed in more detail below but Fig. 2 shows the values of m 
deduced and we see that, provided £o is far enough away from £,™ ax , the values 
obtained from the different lattices are in very good agreement. The maximum 
possible value of £ 0; £o iax f° r behaviour of the form of (8) is obtained by putting 
m = and these values are indicated for the different lattice sizes by arrows 
in fig.l. As £o — » £cT ax > ff^ — ¥ 00 an d resolution for m is lost no matter how 
accurately £o is known; ultimately, to determine m closer to the transition we 
must use larger lattices. The 64 2 data for k < 0.8 can be fitted by the usual form 

m~(K c -K) u (9) 

giving v = .89 ± 0.07 and k c = 0.821 ± 0.005. The fit changes little upon removal 
of points at either end of the range. 

To be confident of the treatment of finite size effects that we have used it 
is desirable to check the assumption (8). Fig.3 shows G» on a 64 2 lattice in the 
region where it is negative together with the best fit of (8) for a number of k values 
close to the transition. The errors on the data are no bigger than the symbol 
sizes but are highly correlated; a correlated error analysis yields \ 2 P er degree of 
freedom of 1.02, 3.18 and 1.75 for k = 0.76, 0.79 and 0.81 respectively. At large 
enough distances, the correlation functions are dominated by the L(p)~ l piece in 
(8) which leads to G\\ ~ — £~ 2 implying |2C], |[ that du = oo. Up to k — 0.81 



our data is in very good agreement with this behaviour of G\\ confirming that 
du = oo for k < k c . ( We would also point out that it is only because of the 
Fourier acceleration in our simulation algorithm that we are able to measure G\\ 
so well at large distances.) The arrow shows £™ ax if the two point function takes 
the form of (8) and we see that the intercept For k = 0.82 the intercept falls well 
beyond £™ ax and (8) no longer accounts for the long distance behaviour of G». 

J It is not practicable to analyse (X(jJ).X(— ft)) directly because it can contain polynomial 
contributions which only affect short distance physics. Taking the Fourier transform and study- 
ing large distances effectively filters out these polynomial contributions. 
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That the two point function suddenly switches from following (8) to something 
completely different is as good evidence as any for the crumpling transition. 

At shorter distances, the action S (1) certainly induces non-trivial interac- 
tions which leads to a discrepancy between the measured G\\ and (8) even in the 
crumpled phase. This is illustrated in fig.4 where, again, the errors are no bigger 
than the symbols; continuing the fits of fig. 3 to shorter distances the \ 2 values 
grow to O(10 3 ) if the very short distance points are included. 

In fig. 5 we plot the specific heat 

C=^((S 2 )-(S) 2 )-l (10) 

for the various lattice sizes. As can be seen the peak height grows steadily, and 
the location of the peak moves toward smaller k, with increasing system size N, 
a classic indication of a second order phase transition. We can try to fit the 
standard divergent behaviour 

C = a + b(K c - K)- a + ... (11) 

to the specific heat for a given system size. On lattices of size 32 2 and smaller 
such a fit cannot be made to our data for any range of k that includes enough 
data points and for N = 32 x 64 we have only obtained points in the immediate 
vicinity of the peak. On the 64 2 lattice a sensible fit can be made to data for 
k < .795 yielding a = 0.17 and k c = 0.819. The specific heat is diverging rather 
slowly and an error analysis for a is complicated by the fact that k c moves quite 
a lot for small changes in a. However, as discussed above, n c is also constrained 
by the mass gap behaviour. A joint fit to the specific heat and the mass gap data 
yields a = 0.2 ± 0.15. Alternatively, the data can be fitted by setting 

a = 2-vd-5 (12) 

and a value for 5 extracted which yields 5 = 0.02 ± 0.3. The central value is 
in much better agreement with the scaling relation a = 2 — vd than previous 
measurements |J although the error shows the need for higher statistics for C . 
Finite size scaling predicts |21| that the maximum value of the specific heat 



behaves asymptotically as 

C = a+ b'L w + ... (13) 

where u) — a/v and L is the linear size of the system. Analysis of our results so far 
is made difficult by the two different shapes of lattice. Applying (13) to the square 
lattice data and assuming the scaling relation a = 2 — ud yields a in the range 0.42 
to 0.66 and v in the range 0.66 to 0.79 with central values a = .56, v = .72. Taking 
L to be the geometric mean of the lattice sides for the asymmetric lattices and 
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fitting all the lattice sizes gives a = 0.47 ± .1, v = .76 ± .05 Q. A previous exercise 
of this kind @ on 16 2 , 24 2 and 32 2 lattices found a = 0.44 ± .05, v = .78 ± .02 
which is in fair agreement. However, these results contradict those obtained from 
the largest lattice alone. Taken together with the absence of a fit of (11) to C 
for the smaller sizes this does suggest that most of the small lattice data is not 
in the asymptotic regime. If we had data for more square lattices this hypothesis 
could be tested by looking for corrections to (12) but the Fourier acceleration 
algorithm limits us to lattice sizes which are powers of 2. 

In terms of the comparison with the KPZ result discussed above we are not yet 
able to draw a definite conclusion but it is intriguing that KPZ is consistent with 
what we know about the critical exponents. If model 1 has specific heat exponent 
oc\ > 0, which the numerical results indicate, then the scaling dimension of the 
energy operator satisfies A e < 1/2 and the KPZ formula with c < 1 predicts 
that for model 2 A f > 1/2 and hence that a<i < 0. This is consistent with the 



most recent simulations of model 2 [15, 16 1, which seem to find that the specific 



heat does not diverge at the critical point, and with the direct determinations of 
central charge for model 1 || |7]]. However, the measurement of c e ff < 1 at the 
critical point is not consistent with proposed modifications of Zamolodchikov's 
c-theorem for non-unitary theories in which the effective number of degrees of 
freedom are supposed to decrease going from the ultra-violet to the infra-red (see 



2^ for a discussion). At k = 0, which is an IR fixed point, we can solve the 
model and hence know that c = c e ff = 3 while at k — k c , which is a UV fixed 
point, c e ff < 1. On the other hand, our simulations yield two point functions 
which indicate that at large distances the crumpled phase of model 1 is described 
by an effective action 

S eff = J~m(K) 2 X- AX + X- A 2 Xd 2 £ (14) 

with m{n) — > at the critical point (it may be that (14) is in fact modified 
by interaction terms and we are at present gathering statistics to try to check 
this). It is straightforward to calculate the correlation function for the energy- 
momentum tensor to find that each embedding dimension contributes 2 to c and 
hence in our case c = c e // = 6 which is certainly consistent with the c-theorem. 
There is more work to be done in understanding the structure of these models 
and their relation to general results in two-dimensional field theory. 



JFW acknowledges valuable conversations with J.Ambj0rn and A.Tsvelik. Some 
of the i860 assembler routines were supplied by C.Michael. This research is 
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Figure Captions 

1. The zeroes of the correlation function G\\ as a function of k for different 
lattice sizes. 

2. The mass gaps deduced from the data in Fig.l by using (8) to analyze the 
finite size effects. 

3. The correlation function G\\ at large distances for different k values on a 
64 2 lattice. The dashed lines show a fit of the form of (8). 

4. The correlation function G\\ at short distances for different k values on a 
64 2 lattice. The dashed lines show a fit of the form of (8). 

5. The specific heat as a function of k for different lattice sizes. The line is 
the best fit of (11) to the 64 2 results. 
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